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The shapes of invariant differential cross section for charged particle production as function of 
transverse momentum measured in heavy- ion collisions are analyzed. The data measured at RHIC 
and LHC are treated as function of energy density according to a recent theoretical approach. The 
Boltzmann-like statistical distribution is extracted from the whole statistical ensemble of produced 
hadrons using the introduced model. Variation of the temperature, characterizing this exponential 
distribution, is studied as function of energy density. 



I. INTRODUCTION 

Inclusive charged particle distributions have been stud- 
ied for a long time to derive the general properties of 
hadronic interactions at high energies. A large body 
of the experimental data on charge particle production 
spectra in baryon-baryon, gamma-baryon and gamma- 
gamma collisions has been accumulated during last forty 
years. However, the underlying dynamics of hadron pro- 
duction in high energy particle interactions is still not 
fully understood. 

Recently, a new unified approach to describe the par- 
ticle production spectra shape was proposed pQ. It was 
suggested to approximate the charged particle spectra 
as function of the particle's transverse momentum by a 
sum of an exponential (Boltzmann-like) and a power law 
statistical distributions: 
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where Erkin = v^r + M 2 — M with M equal to the pro- 
duced hadron mass. A e , A, T e: T, N are the free param- 
eters to be determined by fit to the data. The detailed 
arguments for this particular choice are given in [1 j. For 
the charged hadron spectra a mass of hadrons is assumed 
to be equal to the pion mass. 

Therefore, the hadroproduction process in baryon- 
baryon high energy interactions could be decomposed 
into at least two distinct parts. These parts are char- 
acterized by two different sources of produced hadrons. 
The first one is associated with the baryon valence quarks 
and a quark-gluon cloud coupled to the valence quarks. 
Those partons preexist long time before the interaction 
and could be considered as being a thermalized statistical 
ensemble. When a coherence of these partonic systems is 
destroyed via strong interaction between the two colliding 
baryons these partons hadronize into particles released 
from the collision. The hadrons from this source are 
distributed presumably according to the Boltzmann-like 
exponential statistical distribution in transverse plane 
w.r.t. the interaction axis. The second source of hadrons 



is directly related to the virtual partons exchanged be- 
tween two colliding partonic systems. In QCD this mech- 
anism is described by the BFKL Pomeron exchange. The 
radiated partons from this Pomeron have presumably a 
typical for the pQCD power-law spectrum. Schematically 
figure [I] shows these two sources of particles produced in 
high energy baryonic collisions. This explanation is qual- 
itative, however. This simple model is naive, though it 
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FIG. 1. Two different sources of hadroproduction: red ar- 
rows - particles produced by the preexisted partons, green - 
particles produced via the Pomeron exchange. 



allows to make a number of predictions which have been 
checked experimentally [3H5]- 

In this paper an attempt to study hadroproduction in 
heavy-ion collisions according to the introduced model is 
taken. In these collisions, due to a large number of collid- 
ing partons, extremely high energy densities, comparing 
to those in pp-collisions, can be obtained. The experi- 
mental data measured in AuAu collisions at y/s = 200 
GeV/N by PHENIX [6] and PbPb collisions at = 2.76 
TeV/N by ALICE [7] are considered here. Since the 
centre-of-mass energies per nucleon in these experiments 
are varied by a factor w 14, a unified approach consider- 
ing the energy density is suggested. The energy density 
in heavy ion interactions is known to depend not only on 
the centre-of-mass energy, but also on the centrality of 
the collision. Hence, while the energy densities that can 
be reached at RHIC [6] and at LHC [7 differ significantly, 
the energy density in central collisions at RHIC might be 
of the same order as that in peripheral collisions at LHC. 



2 



Therefore, we use a simple parameterization [8] for 
the initial energy density which is motivated by several 
model calculations. 



e = e (-) a ' 2 N coU P, 
so 



(2) 



where So = lGeV/fm 3 , a « 0.3, /3 « 0.5 and 
200 GeV [8 . Here the second factor is responsible for 
the incident energy dependence, y/s is the cm. collision 
energy, and the third one shows the dependence on the 
number of binary parton-parton collisions N co u which is 
related to the centrality of the collision. 
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FIG. 2. Temperature of charged particles released in heavy- 
ion collisions as function of energy density. 

According to the introduced model, the exponential 
term stands for the radiation of thermalized particles 
with distributions similar to Boltzmann-like thermody- 
namics. Fitting the introduced formula (eq. [TJ to the 
measured experimental data allows to extract only the 



Boltzmann-like statistical distribution from the whole 
statistical ensemble of produced hadrons. This expo- 
nential distribution is characterized by a parameter T e 
analogous to the temperature in classical thermodynam- 
ics, which value is obtained from the fit. Therefore, it is 
interesting to study, how this temperature T e vary with 
the energy density obtained in the collision. 

Figure [2] shows this temperature T e as function of 
energy density. As it was expected, the energy density 
obtained in central collisions at RHIC is similar to those 
in peripherial collisions at LHC, therefore a smooth tran- 
sition between these two experiments is observed. One 
can notice the interesting behavior of the temperature 
as function of energy density (T e oc e 1 / 4 ), which is in 
a good agreement with the Stefan-Boltzmann law. An- 
other observation on the temperature of the final state 
particles is that for high energy densities T e reaches a 
certain limit. This can be explained from QGP theory 
that considers the phase transition temperature T c from 
QGP to hadrons: the system cools until it reaches the 
critical temperature, thus, the temperature of the final 
state should be always below T c . Indeed, for high values 
of e one can notice, that the oberved critical temperature 
is T c « 200 MeV, that is similar to previous theoretical 
observations [9]. 

In conclusion, the experimental results on charged 
hadron production in heavy-ion collisions obtained in the 
PHENIX and ALICE experiments have been analyzed in 
the framework of the new approach. This approach al- 
lows to extract a part of the whole statistical ensemble 
of produced hadrons, those described by the Boltzmann- 
like statistical distribution only. This exponential distri- 
bution is characterized by a parameter T e analogous to 
the temperature in classical thermodynamics. It is found 
that the parameter T e depends universally on the colli- 
sion energy density rather than on the collision energy. 
Finally, the observed behavior of T e is formally similar 
to the Stefan-Boltzmann law well known in classic ther- 
modynamics. 
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